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THE USE OF THEMATIC MAPPER DATA FOR LAND COVER DISCRIMINATION 


PRELIMINARY RESULTS FROM THE UK SATMaP PROGRA^DIE 


1. SUMMARY 

The principal objectives of the UK SAITIaP programme are to determine 
Thematic Mapper (TM) performance with particular reference to 
spatial resolution properties and geometric characteristics of the 
data* Since no data of the UK test sites have been received so far» 
analysis has been restricted to Images from the U.S. and has 
concentrated on spectral and radiometric properties. The results 
indicate that the data are Inherently three dimensional compared 
with the two dimensional character of ^fSS Preliminary classif- 

ication results indicate the importance of the near Infrared band 
(TM 4), at least one middle infrared band (TM S or TM A) and at 
least one of the visible bands (preferably either TM 3 or TM 1). 
The thermal infrared also appears to h<3v.: discriminatory ability 
despite its coarser spatial resolution. For band 4 the forward and 
reverse scans show somewhat different spectral responses in one 
scene but this effect Is absent in the other analysed. From examin- 
ation of the histograms It would appear that the full 8-bit 
quantization is not being effectively utilised for all the bands. 

2. THE SAT!IaP PROPOSAL 

The objective of the SATMaP proposal (Spatial ^alysls of Thematic 
>tapper Products) proposal is to pararaeterise the Thematic Mapper (TM) 
performance with respect to Its spatial resolution properties and its 
geometric characteristics. This assessment will be made in relation 
to the TTf's absolute performance and by comparing results with those 
obtained from the Multlspectral Scanner System (MSS). Three different 
types of measure of spatial resolving capability will be estimated 
these being the modulation transfer function and frequency response 
of the Thematic Mapper, the size of the effective resolution element, 
and the minimum classifiable area for typical land cover types. 

The geometric accuracy tests will be undertaken using large scale 
(1:1250 and 1:2500) Ordnance Survey Maps along with aerial Imagery 
and ground map to spectral band, band to band, and scene to scene 
accuracy. 


3. PROGRESS-TO-DATR 


>htch of tt^ detailed investigation of both the spatial and geometric 
properties of TH data will require imagery of the UK test sites. 
This should become available from NASA following the launch of the 
first Tracking and Data Relay Satellite System (TmSS) and from the 
European Space Agency. Preliminary work has been possible, however, 
based on the availability of 1>l data for Arkansas (Path 23, Row 35, 
8-22-82 ID: 40037-16031) and Detroit (Path 20, Row 31, 7-25-82 ID: 
40009-15A13) and also using data from an 11 channel airborne scanner 
hired from Daedalus. This scanner, which provides the opportunity 
of simulating all seven of the TM channels, was flown over the 
proposed UK test sites in September 1982 and plans exist to coord- 
inate future flights with overpasses of Landsat-4. The analysis has 
concentrated on an examination of the histograms of subscenes, on 
within band variations of spectral response, on interrelationships 
between the bands, on the differing spectral responses of certain 
cover types, and on the separability of these cover types. 

4. DATA PROCESSING 

The analysis of the data was undertaken on I^S System lOl Model 70 
linked to an HP3000 mini-computer. An Optronics C4300 was used for 
direct colour output. The large block size ( 28k bytes}* generated 
by storing three lines of imagery per record is difficult to handle 
on mini-systems and the first stage of processing was to reformat 
the data. Additional basic work included setting up seqv-nces on 
the I^S which displayed the seven channels and selected ti ree band 
combinations in a loop mode. Similar loops were set-up for display- 
ing two dimensional histograms and for moving from viewing the full 
TM image in a subsampled mode to snail areas under the zoom option. 

5. HISTOGRAMS OF SPECTRAL BAJJDS 

Histograms of selected areas for both scenes were examined. Examples 
of these histograms are shown in figure 1 and represent a 512 x 512 
area centred approximately on New Madrid. This area Includes a loop 
of the Mississippi River, a small urban area, agricultural land, and 
woodland and is thus typical of much of the full scene. It is 
notable that the visible bands display a narrow range of digital 
valiies; for example, for band 2 over 957 , of the pixels are found 
within 25 out of 255 digital counts. The histograms are better 
balanced for bands 4, 5 and 6. For sample for band 4 the 95Z limits 
cover 120 digital counts. For classification purposes it is not 
considered that the eight-bit quantization available is being 
effectively utilised even if one takes into account the fact that 
only a single image is being examined. Finally it is of interest to 
note that the data distributed in the near and middle infrared bands 
Is trl-modal but only bl-modal in the visible bands Indicating the 
importance of the former channels for classification purposes. 

*Thls varies slightly between the B, A and P CCT's. 
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6. WITHIN BAND VARIATION IN SPECTRAL RESPONSE 


In assessing che accuracy of classification techniques for Thematic 
Mapper data the consistency of the detector-to-detector response is 
critical. Preliminary studies were undertaken, therefore, to assess 
the significance of this factor for the TM. Ilie approach taken has 
been to examine the spectral response of water bodies where spectral 
frequency responses are generally low and tdiere for wavebands longer 
than 0.7tsB reflectance is both very low and stable. The results 
obtained suggest the existence of anomalies especially in band 4 
for the Detroit scene. Figure 2 shot« a profile of band 4 data over 
Lake Erie. The striping is related to differences in forward and 
backward scans of the sensor. The average difference is approximately 
two digital counts and is clearly present from the eastern edge of 
the Detroit scene to the lake shore near the centre of the image. 
Examination of Reel foot Lake in the eastern portion of the Arkansas 
scene fails to reveal a similar response, nor is it apparent in any 
part of the Mississippi. Either a change in response of the sensor 
or a difference in the data processing is therefore Indicated. 


7. INTER-RELATIONSHIPS BETWEEN BANDS 

The correlations between che thematic bands are shown in Table 1 for 
three sub-scenes. The first two are from the Arkansas scene. The 
first sub-scene is the one described above , and the latter includes 
an agricultural area ami a large region of woodland. The 512 x 512 
sub-scene from Che Detroit scene is centered on Toledo, and includes 
a wide diversity of urban land cover types, agricultural land and a 
portion of Lake Erie. At Che time Che latter was imaged Che cooled 
detectors (for bands 5, 6 and 7) were not operating. For all three 
sub-scenes there are strong correlations between Che three visible 
bands, Che coefficients exceeding 0.9 in all cases. Of the other 
band combinations only bands 4 and 5, 5 and 6, and ft and 7 show any 
strong relationships. Negative correlations are found consistently 
between band 4 in che near infrared and the visible bands 3 though 
as Table 1 shows, che strength of these relationships varies substant- 
ially between the 3ub-scenes. The overall structure of the relation- 
ships can be examined by principal component analysis, the Eigenvalues 
for che three sub-scenes being shown in Table 2. For Che first two 
sub-scenes there are apparently three basic dimensions of variability 
«d)ereas for the third, this is reduced to two dimensions, since bands 
5, 6 and 7 are missing. 


8. CHARACTERISTICS OF INDIVIDUAL LAND COVER CATEGORIES 

In order to examine the utility of Che Thematic Mapper data more 
carefully, six different land cover classes at approximately Anderson 
level 1 were selected frewa the first sub-scene. These included an 
area of water from the sediment-laden Mississippi, woodland, agrlcult- 


5 . 



Fig. 3 Spectral Reaponae Of Selected Cover Type a Prom Sub-scene 1 
Note that band 7 Is the thermal band. 
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Table 1 


CORRELATION ON MATRICES FOR SELECTED LANDSAT 4 THEMATIC MAPPER SUB-SCENES 












1 

2 

3 

4 

5 

6 

1 

1.00 






2 

.967 

1.00 





3 

.953 

.957 

1.00 




4 

-.512 

-.473 

-.631 

1.00 



5 

-.136 

.067 

-.245 

.713 1 

.00 


6 

.244 

.310 

.166 

.362 

.883 

1.00 

7* 

.347 

.391 

.298 

.008 

.567 

.745 

Correlation Matrix 

^qr Sub-scene _2 




1 

2 

3 

4 

5 

6 

1 

1.00 






2 

.950 

1.00 





3 

.935 

.904 

1. 00 




4 

-.148 

-.029 

166 

1.00 



5 

.089 

.207 

.088 

.849 1 

.00 


6 

.353 

.431 

.361 

592 

.904 

1.00 

7* 

.477 

.485 

.455 

.076 

.452 

.674 

Correlation 

_Mat£l_x 

For 

^8_c_ejiq 




1 

2 

3 




1 

1.00 






2 

.923 

1.00 





3 

.989 

.948 

1.00 




4 

-.386 

-.330 

-.210 

Note 

band 

7 la the thennal band* 
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Eigenvalue 


Eigenvalue 


Eigenvalue 


Table 2 


^(WAL 

POOR 


EK7.NVALUES P(» SELECTED SUB-SCENES 


t 2 

.745 .204 

1 2 

.868 .099 

i 2 

.822 .175 


Sub- scene 1 
3 4 

.042 .004 

Sub- scene 2 
3 4 

.025 .003 

Sub- scene 3 
3 4 

.003 .000 


5 6 

.002 .001 

5 6 

.001 .001 


Page 
quality 

7 

.000 

7 

.000 


Note that band 7 is Che thermal band 


ral land and urban land. A "pluiac" class was also selected «dilch 
Includes the pltime of smoke emanating fr<an the power station ami 
drifting over the Mtsslssipi river. The means and standard 
deviations for this class are stwwn in Table 3 at^ the spectral 
response of each class is sumaarised in Figure 3. For the first 
three bands the overall fora of the spectral response is remarkably 
similar, whereas for the near and middle infrared bands considerable 
differences in overall response are found. For the thermal infrared 
band, the means of the classes are closely clustered but It is worth 
noting from Table 3 that the standard deviations are also very low 
suggesting this band may have some discriminatory power. 


9. CORRRUTION STRUCTURE OF INDIVIDUAL COVER CATEGORIES 

Table 4 shows the correlation between the T>] bands for the individual 
cover categories. Considerable differences can be seen between 
their correlation structures ai^ that of the whole sub->scene (Table 
1). For example for water and woodland the correlations between the 
three visible bands are much weaker; the usual weak to moderate 
negative relationship between the visible bands and the near infrared 
band 4 is replaced by strong positive correlations for the imiustrlal 
and plume categories. For the agricultural land and urban land 
categories bands 5 and 6 have moderate to strong relationships with 
the other bands except for band 4, whereas the first two sub-scenes 
as a whole (Table 1} and the water and woodland classes of sub-scene 
1 show much weaker relationships. 

The underlying structure of the relationships is shown in Table 5 
for each of the categories. Except for the water and plume 
categories the eigenvalues are similar to those of the complete 
sub-scenes indicating that there are three basic dimensions. For 
the plume category there are only two dimensions and for the water 
category apparently as many as seven. The latter may simply stem 
from a near spherical distribution of points of this category in the 
residual sub-space. If the eigenvectors of the first two principal 
components are examined (Table 4), it is apparent that although the 
dinenf Ic-nality of the data is similar for the whole sub-scene and 
several of tlie individual land cover categories, the principal 
components are orientated very differently within the seven dimensio- 
nal feature space. In other words different bands contribute to very 
varying degrees to these first two components. Relatively speaking 
bands 2 and 3 contribute least to the first two principal components, 
and bands 4, 1 and S the roost. These results strongly suggest that 
although there may be only three basic dimensions within the data, 
if we tlsh to depict the variability within as well as between broad 
cover categories then more than three spectral bands will have to be 
used. 
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Tiible 3 

MEANS AND STANDARD DEVIATIONS OF COVER CUSSES PROM SUR-SCENE 1 



1 

2 

3 

4 

5 

6 

7 

WATER 

7tean 

87. 4 

38.8 

45.4 

38,0 

7.5 

4.4 

133.0 

Standard 

Deviation 

(1.8) 

(1.0) 

(1.2) 

(2.0) 

(1.4) 

(1.3) 

( .8) 

WOODLAND 

^!ean 

69.7 

26.9 

21.1 

82.4 

51.0 

14.8 

131.5 

Standard 

Deviation 

(2.2) 

(1.1) 

(1.1) 

(5.0) 

(7.6) 

(2.7) 

(1.8) 

AGRICULTURAL LAND 
Mean 79.2 

34.6 

31.6 

112.5 

87.5 

32.1 

140.0 

Standard 

Deviation 

07.9) 

(11.0) 

(14.7) 

(23.7) 

(10.2) 

(9.3) 

(4.0) 

URBA.V 

Mean 

90.1 

39.6 

40.0 

75.5 

76.6 

35.9 

150.0 

Standard 

Deviation 

(11.4) 

(7.1) 

(10.6) 

(12.2) 

(13.1) 

(10.2) 

(4.0) 

PLUME 

Mean 

111.4 

48.9 

55.0 

45.0 

9.5 

5.4 

133.0 

Standard 

Deviation 

(8.5) 

(3.7) 

(3.6) 

(2.2) 

(1.0) 

(1.1) 

(1.4) 

INDUSTRIAL 

Mean 

144.8 

67.3 

77.6 

68.5 

104.9 

67.9 

153.9 

Standard 

Deviation 

(34.0) 

(18.5) 

(23.1) 

(17.2) 

(24.4) 

(18.2) 

(6.6) 


Note that band 7 la the thernal band 
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Table 4 


CORRELATION STRUCTURE OF COVER CATEGORIES 



i 

2 

3 4 5 6 


1 

2 

3 

4 5 6 

2 

.4 

WATER 

2 

.4 


WOODLAND 

3 

.5 

.6 


3 

.5 

.7 



4 

-.4 

-.3 

-.2 

4 

-.1 

.1 

-.1 


5 

.3 

.4 

.4 -.1 

5 

-.2 

.3 

.1 

.6 

6 

.3 

.2 

.3 -.1 .5 

6 

-.2 

.3 

.2 

.5 .9 

7 

.0 

-.1 

-.2 .1 -.2 -.1 

7 

-.2 

.2 

.1 

.5 .A .8 


2 

.9+ 

AGRICULTURAL LAND 

2 

.9+ 

URBAN 



3 

.9+ 

.9+ 



3 

.9+ 

.9+ 



4 

-.2 


-.3 


4 

-.2 

-.1 -.2 



5 

.5 

.5 

.6 

-.1 

5 

.7 

.7 .7 

.3 


6 

.7 

.7 

.8 

-,4 .8 

6 

.8 

.8 .9 

-.2 .8 


7 

.5 

.5 

.5 

■•.7 .3 .6 

7 

.5 

.4 .5 

-.6 .1 

.4 


2 

.9 

PL19E 


2 

.9-»* 


INDUSTRIAL 

3 

.9 

.9 



3 

.9+ 

.9+ 



4 

.9 

.9 

.9 


4 

.8 

.8 



5 

.5 

.5 

.5 

.5 

5 

.7 

.7 

.7 

.7 

6 

.1 

.1 

.1 

.1 .3 

6 

.6 

.6 

.6 

.4 .9 

7 

-.3 

-.4 

-.3 

-.3 -.2 -.1 

7 

-.2. 

- •» 

♦ 

2 

-.1, -.4 -.4 


-t- indlcatee correlation above 0.95 


Table 5 



EICENVALUf;S FOR 

INDIVIDUAL 

LAND 

COVER 

CATEGORIES 



1 

2 

3 

4 

5 

6 

7 

WATER 

.40 

.25 

.13 

.08 

.06 

.04 

.03 

WOODUND 

,77 

.14 

.06 

.01 

.01 

.01 

.00 

ACRI LAND 

.59 

.33 

.06 

.01 

.00 

.00 

.00 

URBAN 

.66 

.25 

.06 

.02 

.01 

.00 

.00 

PLUME 

.94 

.02 

.01 

.01 

.01 

.01 

.00 

INDUSTRIAL 

.82 

.12 

.04 

.01 

.00 

.00 

.00 


Note that band 7 ia the thenaal band- 
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FIRST 2 EIGENVECTOR/RAND LOADINGS FOR SACK COVER CATEGORY 
(Absolute values above O.S are underlined) 

Band number 


Eigenvector 

1 

2 

3 

4 

5 

6 

7 

WATER 1 


.29 

.36 

-.53 

.35 

.25 

-.00 

2 

-.39 

-.04 

-.14 


-.29 

-.21 

.12 

WOODLRNO 1 

.06 

-.03 

-.01 

-.42 

-.85 

-.27 

-.16 

2 

.03 

-.05 

-.05 

.JO 

-.37 

-.15 

-.15 

AGRICULTURAL 







UND 1 

-.53 

-.33 

-.46 

.*J.2 

-.13 

-.26 

-.06 

2 

.3S 

.23 

,25 

.84 

.17 

.05 

-.06 

URBAN 1 

-.4S 

-.30 

-.46 

.03 

-.J2 

-.44 

-.05 

■» 

.U 

.05 

.14 

=,• 3 ?. 

.39 

.05 

.13 

PLl-ME 1 

-.84 

-.36 

-.34 

-.20 

-.05 

-.01 

.00 

*1 

-.20 

.21 

.22 

.06 

.09 

-.04 

7_.J2 

INDUSTRIAL 








1 

-.63 

-.35 

-.43 

-.28 

-.39 

-.26 

,00 

2 

-.35 

-.18 

-.21 

-.12 

.67 

-•17. 

-.05 

WHOLE 1 

-.12 

-.06 

-.14 

.72 

.63 

.21 

-.00 

SUB- 

I7U(7 2 

-.26 

-.2 

-.33 

.50 

-.54 

-.48 

-.12 


Note that band 7 is the thermal bai^. 
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DtSCftIMINATION RF.1VECN COVER CATEGORIES 


The eleeeif tcetory pocenciel of the Themeclc Mapper hands nay be 
exatiined by calculating the divergence bet%reen the claasea^. Table 
7 ahowa the four best bai^a froa the aeven 1>l baMs in dlacrinlnatlng 
between any one class from any of the other classes and the four 
best bands In simultaneously distinguishing between the classes. 
It is notable that except for hand 2 all the bands are It^lcated as 
having significant discriminatory potential for at least two of the 
classes. Overall the results indicate tlM particular significance 
of hands 3, 4 am! S in discrimination between the claSMS. 

Additionally it is interesting that the thermal band, band 7, is of 
independent value in discrimination despite its low spatial resolut** 
ion at^ a less than optimal tine for s«^nsing of this band. 


11. IMPACTS ON UNO COVER CUSSIPICATION 

Whilst the above comments have naturally stressed the anomalies and 
error factors associated with the T>l it should be clearly stated 
that the overall performance is considered to be excellent. Consider- 
able potential exists with the 1>1 data for improved classlficaLion 
performance over that available from 7SS. This is due not only to 
the improved spatial resolvinf capability but to the greater real 
dimensionality of the data. 

The improved land cover classification accuracy will not be achieved, 
however, merely by applying the ti5dltonal per-pixel classifiers so 
frequently used with MSS data. In facr in some circumstances such 
procedures could result in reduced accuracy-. Attention is 
therefore being given to improved techniques and procedure for 
classification of TM data. The Thematic Information Services (TIS) 
of NERC incorporates the Experimental Cartography Unit (ECU) vdiich 
over the last 13 years has acquired large quantities of digital 
vector data and developed s large body of software for manipulating 
and overlaying vector/ polygon data acta. Software haa now been 
written to allow any of theae data contained within the cartographic 
data-base to be Imported and registered with the aatellite dats^. 
In addition the airborne scanner data acquired in September 19R2 is 
aeen as the first acquisition of s regular progrmma of flying and 
again the imagery is being input si^ registered with the map and 
satellite data. The opportunity therefore exists to apply 
classification techniquaa which exploit per-pixal, textural and 
contextual algorithms within tlM framework of a probabilistic tree 
classifier and an integrated multi-level and multi- pa rams tar data 
set. As for the characterleatlon of the TN geometry and spatial 
rasolvit^ propartiaa tlgnlf.cant further prepress aiMits TM data of 
the UK but preliminary results uaing the airborne Daadalaa and MSS 
data are promising and arc being applied in a number of etudiee 
ranging fr<Mi ecology to epidemiology. 
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Table 7 


DIVERGENCE BET^TREN CLASSES: BEST FOUR BANDS FOR 
DISCRIMINATING BETWEEN CLASSES 

1 2 3 4 5 6 7 


Water 

Agricultural X 

Land 

Woodland X 

Ur band 

Plune 

Industry 

All classes 


X X X X 

XXX 

X 
X 
X 
X 


XXX 
XXX 
XXX 
X XX 

XXX 


Note that band 7 Is the thermal band 


12. PRELIMXHARY ANALYSIS OP THE FIRST UK S(XNE 


lh€ first TM insgs of the UK mss obtained In January 19S3. This 
image «ms received via X-bat^ transmission at Fucino, Italy and 
processed at ESRIN. Frascati. Unfortunately the scene includes 
nelthar of our UK test areas but does Include tlM towns within 
%ihlch our two research institutes are foimd. At present only a 
single. TM band 4 image has been received but sone preliminary 
evaluation is possible. Two obvious defects are found in the image, 
^le first relates to shifts due to a lack, of registration beti#een 
forward and reverse scans. The strength of the shift is up eight 
pixels though normally it is three or less. This shift appsrently 
arises because precise orbital data have not been used in image 
processing so that corrections for minor changes in satellite 
orientation have not been Included. A number of drop~outs can be 
seen thoughout the image. Their origin is not currently understood. 
Currently steps are heli^ taken to correct these minor defects. 
These have digital counts between 50 and 80 ratter than zero. It 
should be emphasised, however, that ESA released this scene to 
European PI*s at the earliest possible time with the knowledge that 
these problms existed. In terms of Information content this image 
represents s very substantial improvement on landsat MSS data in 
terras of the detectability of ground targets especially in rural 
areas. In the upper and lower parts of the scene, it is apparent 
that individual fields are readily detectable, see figure 4 which 
certainly was not true in this area with MSS data. In several places 
narrow hedgero«ra can be seen which are less than 6in across and most 
of which are less than 3ra across. 

The River Thanes can be seen meandering across the northern half of 
the image. The large black areas are gravel pits, now filled with 
water. The nain through *’oute in the area (the M4) can be seen in 
the southern part of the . age- Whereas the image is significantly 
better in rural areas than the MSS, in the urban areas the improve- 
ments are much less narked. The urban area is represented by the 
mottled dark tones across the full width of the central part of the 
image. Compared with images of urban areas in the US such as Detroit 
and Washington D.C., the results are apparently rather disappointing. 
However it must be noted that the density of urban development is in 
general much higher in English than American cities and hence much 
less ground detail will be detectable. For resolution of the Jetails 
of urban structure of cities like Reading, spatial resolutions of 
20n or better will be required. 
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